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Abstract 

A model for the condensation of a dusty plasma is constructed by considering the spherical 
shielding layers surrounding a dust grain test particle. The collisionless region less than a collision 
mean free path from the test particle is shown to separate into three concentric layers, each having 
distinct physics. The method of matched asymptotic expansions is invoked at the interfaces between 
these layers and provides equations which determine the radii of the interfaces. Despite being much 
smaller than the Wigner-Seitz radius, the dust Debye length is found to be physically significant 
because it gives the scale length of a precipitous cut-off of the shielded electrostatic potential at 
the interface between the second and third layers. Condensation is predicted to occur when the 
ratio of this cut-off radius to the Wigner-Seitz radius exceeds unity and this prediction is shown 
to be in good agreement with experiments. 
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I. INTRODUCTION 



Condensation of a dusty plasma [1] into a crystalline state was proposed by Ikezi in 1986 
[2] and demonstrated experimentally eight years later by a number of research groups [3-6] . 
The subject has been reviewed by Morfill et al. [7] and most recently an experiment to test 
dusty plasma physics has been set up on the International Space Station [8]. 

The original model [2] for this process was motivated by Monte-Carlo calculations [9] 
which predicted that a Coulomb crystal would form when the Coulomb interaction energy 
between two adjacent charged particles in a one-component-plasma exceeded their thermal 
energy by some factor. The Coulomb interaction energy for charged particles with density 
n is the electrostatic energy of one particle in the potential of an adjacent particle located 
at the Wigner-Seitz interparticle separation distance [10] 
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According to the Monte-Carlo calculations, condensation of charged particles into a crystal 
should occur when 

Z 2 e 2 

r = A * £ 170 ( 2 ) 

where Z is the charge on each particle and T is the temperature of the particles. As noted 
by Ikezi [2], Eq.(2) could be a very poor estimate for dusty plasmas (which are a three- 
component-plasma), but lacking a better model, Eq.(2) has often been used as a benchmark 
for dusty plasma crystallization experiments. The experiments show [6] that the actual 
value of T required for condensation is two to three orders of magnitude larger than that 
predicted by Eq.(2). Thus, while Ikezi's original postulate that dusty plasmas can condense 
into crystals has been experimentally validated, there has not been a quantitative model 
predicting the value of Y necessary for condensation to occur. 

Interactions between adjacent particles in a plasma are intimately related to the concept 
of Debye shielding. According to this concept, any plasma particle can be considered to be a 
test particle surrounded by a screening cloud of adjacent particles. The screening completely 
cancels the test particle field at distances much greater than the Debye length. Screening 
may be accomplished by adjacent particles of either the same polarity as the test charge or 
opposite polarity, but is subject to the constraint that the test particle cannot be moving 
faster than the thermal velocity of the shielding particles [11]. For example, if the test 



particle is an electron, it is shielded by the repulsion of other electrons in the presence of a 
uniform neutralizing ion background, but it is not shielded by ions because it is moving too 
fast for ions to respond. On the other hand, ions are shielded by both electrons and ions. 
In a dusty plasma one might thus reasonably expect dust grains to be shielded by electrons, 
ions, and other dust grains. 

The standard model of Debye shielding is based on the Boltzmann relation, an equilibrium 
solution to the fluid equation of motion for each species a such that the force due to the 
electric field balances the force due to the gradient of an isotropic scalar pressure, i.e., 

= -n CT g CT V0 - VP CT . (3) 

Three critical assumptions are intrinsic to the standard model of Debye shielding, namely: 
(i) it is assumed that the plasma is sufficiently collisional that the concept of an isotropic 
scalar pressure P a = n a KT a is valid, (ii) it is assumed that a Boltzmann dependence n a = 
n a oexp(—q a (f)/ nT a ) exists relating the local density n a to the system-averaged density n a0 , 
and (iii) it is assumed that \q a (f> / nT a \ « 1 so that the Boltzmann relationship may be 
linearized giving n a /n a Q = 1 — q a (f)/ KT a . The standard model for Debye shielding of a test 
particle with charge qr results when the linearized Boltzmann relationships of the various 
species are substituted into Poisson's equation giving the Yukawa-type solution 0(r) = 
§Texp(— r / \d) / Ane^r where 



and 
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AL = • (5) 

The summation in Eq.(4) is restricted to species that participate in the shielding and so 
excludes all species having thermal velocity slower than the species of the test particle. 

When r — > 0, the Yukawa solution diverges, violating assumption (iii) that \q a <j)/ KT a \ « 
1 and causing the standard model to fail to be internally self-consistent. This failure of the 
standard model of Debye shielding has been noted previously, e.g., see Lampe, Joyce and 
Ganguli [12] for criticism regarding assumptions (i)-(iii). In addition, Hansen and Fajans 
[13] have shown that trapping can affect Debye shielding in a pure electron plasma while 
Goree [14], Zobnin et al. [15], and Lampe et al. [16] have shown that trapping of ions can 
affect shielding of a dust grain. 
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The issue of how to treat Debye shielding when \q a (j) / KT a \ » 1 is especially critical for 
the dust condensation problem, because \q a (f) / K,T a \ is essentially the same as T. Consensus 
does not exist on how to address this issue. 

Furthermore, the form of Eq. (4) is such that the sum on the right hand side is dominated 
by the term having the smallest X 2 Da and, since dust particles are both cold and highly 
charged, the dust Debye length is typically much smaller than both the electron and ion 
Debye lengths. One might expect that the system Debye length X D should be very nearly 
the dust Debye length A^, but this point of view has usually been rejected. The dust Debye 
length Xod is typically so small that it is less than a, and questions have been raised as to 
whether such a short shielding length has physical significance since the standard Debye 
argument is based on the implicit assumption that there is a statistically large number of 
particles in a sphere having the Debye radius. This is clearly not true for dust particles 
in a sphere with radius X^d if X^d is less than a. Nevertheless, Wang and Bhattarcharjee 
[17] and also Otani and Bhattarcharjee [18] argued that some sort of shielding does occur 
at the scale of Xod but the only support for this point of view was demonstration [18] of 
some time-averaged correlation effects at the scale of Xod in a one-dimensional numerical 
simulation that would only crystallize if artificially annealed. Most other authors ignore 
dust self-shielding on the presumption that the Debye shielding concept does not make 
sense when a Debye length is smaller than a. 

We present here a model for a dusty plasma on the verge of condensation. This model 
takes into account both collisional and collisionless behavior in three dimensional geometry, 
avoids inappropriate use of fluid theory, shows that the dust Debye length has important 
physical significance even though it is much smaller than a, and predicts a condensation 
threshold in good agreement with experimental measurements. The derivation identifies four 
physically distinct concentric regions surrounding a test charge. The ions are collisionless 
in the innermost three regions but collisional in the outer, fourth region (the collisionless 
nature of ions in the inner three regions is consistent with the assumptions inherent in dust 
grain charging theory). The method of matched asymptotic expansions is used to locate 
the two interfaces between the first three regions and knowledge of these interface locations 
is then used to give the criterion for condensation. The T e >> Tj,T^ temperature regime 
of typical dusty plasma condensation experiments is assumed and ions are assumed to 
be singly charged (the theory could be extended to arbitrary temperatures without great 
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difficulty, but this would unnecessarily complicate the model). 

The paper is organized as follows: Section II reviews relevant aspects of dust charging 
theory and sets up a dimensionless parameter space suitable for comparing the model to 
experiments. Section III uses collisionless Vlasov theory to calculate the ion, electron, and 
dust grain densities and shows that when |e0//tTj| > 1, the ion density differs from the 
Boltzmann model; this difference demonstrates the inappropriateness of fluid models in 
this regime and resolves the paradox associated with divergence of the Yukawa solution 
at small r. Section IV shows that the vicinity of a test particle can be divided into three 
concentric spherical regions each having distinct physics determined by the magnitude of 
|e0//cTj|. Section V derives approximate solutions to the Vlasov- Poisson system for these 
three regions and Section VI derives matching conditions across the two interfaces between 
the three regions. Section VII uses the matching conditions to deduce a condition for dust 
condensation and compares the model predictions with experiments. Section VIII provides 
a summary and discussion. 

II. DUST CHARGING AND DUSTY PLASMA PARAMETER SPACE 

Two independent parameters characterize the dust grains in a dusty plasma: the grain 
radius and the Wigner-Seitz radius a. In order to develop a model based on dimensionless 
parameters, the ion Debye length 



will be used as the 'yardstick' by which all lengths are measured. A bar will be used to 
denote lengths normalized to the ion Debye length so that the normalized Wigner-Seitz 
radius, for example, is 



The two quantities a and fa constitute the coordinates for a dimensionless dusty plasma 
parameter space. 

In order to avoid confusing minus signs, the electrostatic potential will be replaced by 
the positive dimensionless variable 




(6) 
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and ipd will denote the potential on the surface of a dust grain. Thus, positive ip attracts 
ions but repels both electrons and dust grains. 

When dust grains are placed in an electron-ion plasma, some fraction of the electrons 
attach to the dust grain surface, causing the dust grains to become negatively charged and 
reducing the density of free electrons. The quantitative theory of dust charging, summarized 
in Ref. [1], combines collisionless Vlasov theory with an analysis of trajectories of individual 
particles as they approach a finite radius charged sphere. The particle trajectories are 
assumed to be governed by orbital-motion-limited (OML) theory [19-21] wherein particle 
trajectories are assumed to be collisionless and completely determined by considerations 
of conservation of angular momentum and conservation of energy. There has been some 
question [6] about the extent to which the standard dust charging model applies to dust 
grains in an electrode sheath, the typical situation for terrestrial dusty plasma condensation 
experiments, but not for the zero-gravity dusty plasma condensation experiment on the 
International Space Station. We assume in this paper that the standard dust charging 
model is applicable so that the effect, if any, of electrode sheaths on dust charging is small. 

The standard dust charging model shows that dust grain charging is governed by the 
dimensionless parameter 

P = Aim d0 \ 2 Di r d = 4irn d0 \ 3 Di f d = — (9) 

(X 

where P has the functional dependence 



P = T ~ ( 1 + T J J^expiW/Te). (10) 
Vd V WdJ V miT e 

Global quasineutrality gives 

Z d n d0 + n e0 = n i0 (11) 

where Z d is the number of electrons captured by a dust grain. We define the electron capture 
factor 

a = Z d n d0 /n i0 (12) 

so that a = 1 corresponds to having all the electrons attached to the dust grains while 
a = corresponds to having no electrons attached to the dust grains. The quasineutrality 
condition, Eq. (11), can thus be expressed as 

a + ^ = l (13) 
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and dust charging theory [22] shows that 

a = Pi> d . (14) 

Since ip d and a are functions of P, they have functional dependence a = a{a,f d ) and 
= ip d (a,f d ). 

Combining Eqs. (9), (12) and (14) shows that 

~ A d ~TT = f d*Pd (15) 

so that Z d becomes large if ip d is finite, f d is not infinitesimal, and Ann^X^ is large. Equation 
(15) is just the normalized version of the potential <pd = —Zde/Aneord of a sphere of radius 
Td with surface charge —Z d e. This result is actually slightly incorrect for a shielded dust 
grain, because, as shown in the next paragraph, the shielding cloud surrounding a dust grain 
depresses the potential at the grain surface to a value below the value given by Eq.(15). 

To understand this potential depression effect due to shielding, consider the potential (f> 
on the surface of a sphere with charge Q and radius r sp here surrounded by a shell of shielding 
charge — Q at radius r s heii- The potential on the surface of the shielded sphere is given by 

^>{r sphere) = / dr— = , (16) 

Joo Or A-KSq \r sphere r sh ellj 

a result obtained by taking into account the contributions to d<f>/ dr from both the sphere and 
its shielding charge. The ratio of the surface potential of the shielded sphere to the potential 
of an identical unshielded sphere is <j)(r S phere)/4>vac — 1 — 'sphere/ f 'shell where <p vac is the surface 
potential of the unshielded sphere. If r s heii ~ r sphere — where is the nominal Debye 
length, then 4>{j- sphere ) / (p vac ~ \ D / [r sphere + -Vd) so the potential of the shielded sphere will 
be greatly depressed from its vacuum value if r sp here » Ad- This indicates that being 
highly charged is insufficient for a dust grain to have a large potential; it also needs to 
have fd « 1. The model of dust charging thus has the implicit assumption that f d is small 
compared to unity and this assumption will be made in the remainder of this paper. 

Figures l(a)-(c) show contours of constant ipd, and Z d / 4:im i0 \ 3 Di as determined by 
Eqs. (9), (10), (14) and (15) for the parameters of Ref. [3], a typical dust crystallization 
experiment. Since dust grains in a given experiment have a fixed ratio f d /a, a specific 
experiment is characterized by a sloping straight line in a, f d parameter space. Moving 
up and to the right along such a line corresponds to making A^ smaller whereas moving 



down and to the left corresponds to making Xdi larger. Densities in an experiment are 
typically measured by Langmuir probes which have an uncertainty of -50%, +100%, so 
that the density of = 10 9 cm" 3 reported in [3] would actually be in the range 5 x 10 8 
cm -3 < rti < 2x 10 9 cm -3 . This factor of four range of densities [4] corresponds to the straight 
line segment labeled 'expt' in Figs.l(a)-(c). This line has a slope given by f^ja = r&ja. The 
left end of this line is the point in parameter space calculated using the lower estimate for 
the density, while the right end corresponds to using the upper estimate for the density. The 
length of this line effectively represents the density measurement error bar. To the extent 
that charging theory is correct, the range of possible values of a, ipd, and Zd/ ^nn^X^ for 
the experiment are given by the intersection of the contours with this 'expt' line. 




4 
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FIG. 1: Dusty plasma parameter space for Chu and I experiment (density range indicated by short 
line labeled expt): horizontal axis is a = a /Xdi, vertical axis is fd = rd/XdV, (a) shows contours of 
constant tpd, (b) shows contours of constant a, (c) shows contours of constant Zd/4:imioX di . 
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III. VLASOV MODEL OF CHARGED PARTICLE DENSITY IN THE PRESENCE 
OF A POTENTIAL 



A typical dust grain will be considered as a test particle inserted in a plasma consisting 
of electrons, ions, and other dust grains. The origin of a spherical coordinate system will be 
defined to be at the center of this test particle. Typical dust condensation experiments have 
neutral pressures ~ 10 2 Pa (corresponding to a neutral density n n ~ 3 x 10 16 cm -3 ). Since 
neutral cross-sections are a ~ 3 x 1CT 16 cm 2 , the mean free path for ion-neutral collisions 
is l rn fp = (n„(x) -1 ~ 1 mm which is at least an order of magnitude larger than shielding 
scale lengths. The last collision experienced by an ion in the vicinity of the test particle 
will have occurred outside a sphere having a normalized diameter of the order of l m f p ] such 
a sphere is shown schematically in Fig. 2 and lies at the interface between regions 3 and 
4. Thus ions can be considered collisionless inside regions 1 to 3 and collisional in region 4; 
the details of regions 1 to 3 will be discussed later. This separation of space into collisional 
and collisionless regions is similar to the arguments used in dust grain charging theory (the 
OML assumption underlying dust grain charging theory is based on angular momentum 
conservation which can only be true if a particle has no collisions). 

Electrostatic potential is undefined with respect to a constant; following convention, 
we choose this constant such that <fi = at infinity. Collisions make the distribution 
function Maxwellian in region 4 and this provides a boundary condition for the collisionless 
distribution function in regions 1 to 3. The distribution function in regions 1 to 3 must 
satisfy the collisionless Vlasov equation and so must be must be a function of constants 
of the motion [23, 24]. The relevant constant of the motion here is the particle energy 
W = m a .v 2 /2 + q a (f), and so the distribution function in regions 1 to 3 is 

, ( v / m a \ 3/2 / m a v 2 /2 + g CT 0(r) \ 

f ^ v) = n - ) exp { ^ — ) ■ (17) 

This is the right choice because / CT (r, v) is not only a function of a constant of the motion 
W but also joins smoothly to the region 4 Maxwellian solution where = and the plasma 
is collisional. 

Electrons and dust grains experience a repulsive force upon approaching the negatively 
charged dust grain test particle and so are slowed down with some particles being slowed 
down to zero velocity and reflecting. Thus, electrons or dust grains near the dust grain 
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region 4 




collisional 
outside 
"this circle 

FIG. 2: Sketch of concentric regions surrounding a dust grain test particle of normalized radius 
r d . Diameter of outermost dashed circle is l m f p , the normalized mean free path for collisions, so 
plasma is collisional in region 4 outside this circle. Regions 1, 2, and 3 are collisionless and have 
interfaces on the dashed circles having normalized radii fj and f a . 



test particle can have zero velocity. The respective electron and dust grain densities in the 
vicinity of the dust grain test particle are thus given by 

f°° ( e6(r)\ 

n e = J 4:7rv 2 dvf e (r, v) = n e0 exp 11 , (18) 

n d = J 4irv 2 dvf d (v, v) = n d0 exp (^-^-^] . (19) 



Using Eq.(8), the electron and dust grain densities normalized to their average values are 



■ exp(-^T t /T e ) (20) 

n e0 



and 



m = exp {-Z^Ti/Ta) . (21) 

These densities are identical to the Boltzmann relation and so demonstrate that collisionless 
kinetic theory agrees with fluid theory for negatively charged particles near a negatively 
charged test particle. 
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Ion behavior is fundamentally different because ions, being positive, are accelerated as 
they approach the negatively charged dust grain test particle; this means that there are no 
zero velocity ions near a dust grain. The slowest ion is one that has fallen into the negative 
<f> well with zero initial velocity at the edge of the well and the velocity of such an ion will 
satisfy 

miV 2 /2 + e<p = 0. (22) 
Using Eq.(8) it is seen that this minimum possible ion velocity can be expressed as 



*W = \ /^V- (23) 

V m 

Evaluation of the ion density in the vicinity of the dust grain test particle therefore 
requires invoking a lower limit at i> m i n for the velocity integration over the distribution 
function. The resulting ion density is thus 



POO 

/ 4irv 2 dvfi(r, v) 



n i0 



Hi 

( -^-Y^exp ($) I" 47rv 2 dvexp ( -^-) . (24) 

By defining £ = v/ ■\f2KT i /m i , the normalized ion density can be expressed as 

ni_ 4e^ 



poo 

/ mexp(-e 2 ) 



= (l-erfV?) + 4^ (25) 



eiiz= 4= / e-« 2 ^ (26) 



where 

erf z = — 
V 71 " 

is the error function. The second line in Eq.(25) is obtained using the identity 

poo 7 poo poo 

/ dt-fcexp {-?)) = ^exp(-e 2 )-2/ « 2 exp(-£ 2 ). (27) 

For small arguments, the error function may be approximated 

2 / r 3 \ 

lim erf z = -r= I z . (28) 



z-»0 7T 

Thus for ^ << 1 and hence << 1, the normalized ion density has the form 

~ 1 + V (29) 
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which is the same as the Boltzmann result given by fluid theory. 
However, because 

lim [l - erf \ = 0, 

the normalized ion density when ip » 1 is 

--^v^; (so) 

nio V 71 " 

this is much smaller than the fluid theory Boltzmann relation prediction that ni/n i0 = 
exp(^). Equation (30) thus demonstrates a failure of fluid theory and its associated Boltz- 
mann relationship when ip » 1. This failure occurs because the concept of ion pressure 
no longer makes sense when ip » 1. The pressure concept is based on the assumption 
that particles have an isotropic Gaussian distribution of random velocities about some mean 
velocity whereas when ip » 1, ions in reality are falling into a deep potential well and 
do not have a random distribution of velocities about some mean velocity. Equation (25) 
and the distinction between its small and large ip limits have been previously discussed by 
Laframboise and Parker [25] in the context of electrostatic probes and by Lampe, Joyce and 
Ganguli [12] in the context of dusty plasmas. 

The lower limit of the integral in Eq.(24) causes the ion distribution to have an r — v 
phase-space 'hole' in the vicinity of the dust grain since f(r,v) = for velocities below 
f m in- It has been argued by Bernstein and Rabinowitz [26], Laframboise and Parker [25], 
and Lampe [27] that for a certain class of radial potential profiles, another sort of phase- 
space hole can also exist. This additional hole results from a rather subtle barrier that 
can occur because for a certain range of the angular momentum J the effective potential 
U e ff(r) = q<f>(r) + J 2 /2mr 2 can have a small local maximum. This barrier prevents access 
to small r by particles having a certain range of W and J. If such a barrier exists, the radial 
ion density profile will differ somewhat from the predictions of OML theory, because the 
ions that cannot pass by this barrier will have a radial turning point at a larger radius than 
predicted by OML. However, since Poisson's equation shows that ip is essentially a double 
integral of the net charge density up to a radius r, changes in the turning point of small 
classes of ions should not have a major effect on the ip profile, i.e., small corrections to the 
OML model should not result in a significant collective effect. Lampe [27] has shown that 
the error introduced by omission of consideration of these centrifugal force barriers is very 
small for dusty plasmas and so we will ignore this correction to OML theory. 
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Another correction to OML theory results from consideration of ion capture by the dust 
grain which also causes a hole in phase-space [20, 26]. As shown in Ref.[26] capture of ions 
by the dust grain reduces the number of ions moving radially outward from the dust grain in 
comparison to the limiting situation where the dust grain does not capture any ions so that 
all ions are perfectly reflected from the dust grain. Taking into account the reduction in the 
number of outward moving ions compared to inward moving ions would require replacing 
the distribution function prescribed by Eq.(17) by a distribution function of the form [26] 
/ = f + + /_ where f+(W) is the phase space density of ions moving radially outwards from 
the dust grain and f-(W) is the phase space density of ions moving radially inwards. If ions 
are perfectly reflected at the dust grain then /+ = /_ in which case Eq.(17) is appropriate, 
but if ions are captured by the dust grain then /+ < /_ and a more complicated prescription 
than Eq.(17) would have to be used. We will assume that the fraction of ions incident at 
f Q which are captured by the dust grain is so small that Eq.(17) is a reasonably accurate 
prescription for the ion phase space density. We are thus assuming that the r projection 
of ion motion has a reflecting trajectory so that there are equal numbers of ions moving 
radially inwards and outwards in the dust grain shielding cloud. This assumption will be 
validated later. 

Finally, we will also ignore ion trapping [14-16] , but will later make some brief comments 
about the extent to which trapping might be important. 

IV. IDENTIFICATION OF THREE REGIONS FOR THE POTENTIAL 

Poisson's equation 



relates the densities of the various species to the electrostatic potential. Assuming spherical 
symmetry about the dust grain test particle and using Eq.(8), Poisson's equation can be 
recast as 




(31) 




(32) 
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Using Eqs.(20), (21), (25) for the normalized densities and also Eq.(12), the normalized 
Vlasov/Poisson system becomes 

hlf ( f2 t) = i 1 - CTf (^)) + ^ - (1 - exp (-f ) - aeM-Z^ 

V v ' v v ' dust 

ions electrons 

(33) 

where Z = Z d Ti/T d is presumed to be large compared to unity since the dust grain is highly 
charged and T; L > T d . Equation (33) is a nonlinear ordinary differential equation for ip and 
is consistent with the collisionless Vlasov equation. Since the densities were obtained using 
the collisionless Vlasov equation, this system will be called the Vlasov/Poisson system to 
distinguish it from the fTuid/Poisson system. 

We now argue that three distinct regions exist for ip such that in each region the 
Vlasov/Poisson system has a different form. The location of these regions is sketched in 
Fig. 2 and, going outwards from the surface of the dust grain test particle, these regions and 
their interfaces are defined by: 

Region 1 is where ip d > ip > 1 and exists because the grain potential ip d is large compared 
to unity [see Fig. 1(a)]. Region 1 is a sheath-like inner region where the ion density has the 
non-Boltzmann behavior given by Eq.(30). 

Region 2 is where 1 > ip > 1/Z d and is depleted of dust grains. Shielding in this region 
is provided mainly by ions. 

Region 3 is where 1/Z d > ip and this region extends to region 4 where collisions set in 
and where the potential goes to zero. Shielding in region 3 is done mainly by dust grains 
and this shielding takes place over a very short characteristic length, causing an extremely 
sharp cut-off of the potential. 

The radii of the respective interfaces between regions 1 and 2 and between regions 2 and 
3 will be called fj and f as indicated in Fig. 2 (the subscripts i and o stand for inner and 
outer interfaces). The values of fj and f will be unknowns to be solved for; determining 
these radii is the crux of the problem. 
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V. APPROXIMATE SOLUTIONS TO THE VLASOV/POISSON SYSTEM FOR 
THE THREE COLLISIONLESS REGIONS 



The three collisionless regions will now be discussed going from the outermost (region 3) 
to the innermost (region 1). 

A. Region 3 solution: ijj < 1/Z 

In region 3, Eq.(33) can be approximated as 

lons v v ' dust 

electrons 

= (l + (l-a)?± + aZjil>. (34) 

Using T e >> Tj, this has the Yukawa- type solution 

^ 3 = jz exp (- VaZ + 1 (f - f )) . (35) 

The coefficient in Eq.(35) has been chosen so that ^ 3 = 1/Z at f = f a . The effective 
shielding length in region 3 is the dust Debye length 

\ Dd = \ Dl /Vl + aZ (36) 

which is much smaller than the ion Debye length since a is of order unity [see Fig. 1(b)] and 
Z = Z d T t /T d » 1. 



B. Region 2 solution: l>ip>l/Z 

The exp (—Ztp) term is dropped from Eq.(33) in region 2 because Zip is large. Taking 
into account T; << T e and ip < 1, Eq.(33) reduces to 

ions electrons 

or 

f 2 df \ df ) ^ a (^) 
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By considering ip + a as the unknown, it is seen that ip + a has solutions of the form 
f _1 exp(±f); the exponentially growing solution is allowed here because region 2 does not 
extend to infinity. A convenient way of expressing the general solution is 

A cosh (f — f ) + B sinh (f — f ) 

ip + a — = . (39) 

r 

In order to have tp = 1/Z when f = r , we choose 

(40) 



and leave B undetermined. Thus 

f Q ( ■= + a J cosh (f — f ) + 5 sinh (f — f ) — ar 

^ 2 = L = (41) 

r 

is the region 2 solution with coefficients arranged so that ^2 — 1/Z when f = f . 



C. Region 1: Inner region 

In this region ip > 1 and we assume that tpTi/T e « 1 so that Eq.(33) reduces to 
1 d 





(42) 



where the electron term has been dropped because 1 — a is significantly less than unity and 
ip is assumed to be larger than unity. Equation (42) can be written as 

F + r¥ = M# (43) 

where 

= -E= « 1. (44) 

Since /i << 1, the right hand side of Eq.(43) may be neglected compared to either of the left 
hand terms in which case the approximate solution to Eq.(43) is the vacuum-like solution 

if) = — = — (45) 
r 

where c and d are constants to be determined. The coefficient d provides for the slight 
depression of the grain potential due to the shielding cloud. The d term is allowed because 
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region 1 is of finite extent and so finite d is not inconsistent with ip vanishing at infinity 
since infinity is not located in region 1 . 

From Gauss' law, the radial electric field E r at the dust grain surface is 

Ane r 2 d E r = -Z d e. (46) 

Since E r = —d(p/dr = (kT^/ eXoi) dip jdf , the boundary condition at the grain surface 
can be expressed as 



dr J - d 4irn i0 \ 3 Di r d 



This gives 

and so, using Eqs.(7) and (12), 



4irn i0 \ 3 Di 



(48) 



aa 3 



c = — . (49) 



By assumption — 1 at r j and so, using Eq.(45), 



in which case 



r 

Thus the region 1 potential is 



%r + d = 1 (50) 



d=l--. (51) 



aa 3 ( _ 1 aa 3 



+ 1 - —-r- f 



^ = V _ - - 3 7 ; (52) 

r 

this satisfies Gauss's law at the dust grain surface and also gives ip — 1 at f = fj. The 
potential on the grain surface is 

aa 3 ( 1 1 \ , . „. 

^ = ^k-7,) + L < 53) 

VI. MATCHING THE SOLUTIONS 
A. Matching principle 



Matching consists of arranging for equality of ip and tp' at the two interfaces between the 
three collisionless regions; the necessity for continuity of ip and ip' across an interface is 
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established by integrating Eq.(33) twice across the interface. Solutions on the left and right 
hand sides of a matching radius f m are of the general form ipieft = L(f)/f and ipright = R(r) / f 
and so matching requires 

L{f m )/f m = 4> m = R{f m )/f m (54) 
L'(f m ) = R'(f m ) (55) 

where ip m = 1/Z when f m = f Q and ip m = 1 when f m = f\. Here L and R are the left 
and right hand numerators: L is the numerator of ipi and R is the numerator of i/j 2 when 
r m = ff, L is the numerator of ip2 and R is the numerator of ^3 when f m = f Q . 



B. Matching of ip2 and ^3 at r Q 

The ^ 2 and ^3 solutions have already been arranged to satisfy Eq.(54) (i.e., ip 2 = ip3 = 
1/Z at f ). The derivative matching condition, Eq.(55), is satisfied if 

B = a-^VaZ +1 (56) 

and so 

f Q ( — + a J cosh (f — f D ) + fa — -^-VaZ + 1 ) sinh (f — f ) — af 

^ 2 = ^ — (57) 

r 

smoothly matches to ip 3 at f D . The actual value of f is undetermined at this stage and will 
be found later. 



C. Matching of ipi and ip2 at rj 

Since Eq.(54) requires R = f\ in order to have ipifi) = 1, Eq.(57) provides the relation 

r» = f + a^j cosh (f, - f ) + (a - aZ + 1 j sinh (f* - f ) - or*. (58) 

The condition L — fi when ■?/> = 1 has already been arranged by the form of Eq.(52). 

The condition V = R', found by taking derivatives of the numerators of Eqs.(52) and 
(57), is 
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1 ~~ = r °\~z + j sinl1 ~~ r °) + ( a ~~ ~z~ v /(yZ + 1 ) cosl1 ^ ~ ~~ a ' ^ 

Equations (58) and (59) constitute two coupled equations in the unknowns fj and f G . 
Using Z >> 1, these equations reduce to 

(1 + a) fi = af Q cosh (f j — f G ) + a sinh (fj — f G ) (60) 

1 + a — = af sinh (f j — f G ) + a cosh (f j — f G ) . (61) 

3rj 

For given a and a these nonlinear equations can be solved numerically for fj and f a . Since 
a = a(a, fd), this means that for any point in a, f^ parameter space, one can calculate a 
and then calculate fj and f Q . Thus, we can consider fj = fi(a,f d ) and r Q = f {a,f d ). 

Once fj and f D are known, the solutions ipi,ip2, and ^3 are all determined and match 
smoothly across the interfaces. The Vlasov/Poisson equation is thus solved all the way from 
the grain surface to infinity. The potential falls off abruptly at f > f with a scale length 
given by the dust Debye length. The dust Debye length is thus of physical importance 
even though it is much smaller than the inter-particle spacing. No paradoxes occur due to 
this situation because the solution for ip is multi-scale and more complicated than a simple 
Yukawa type potential. In particular, the dust shielding does not take place in a sphere 
having a radius equal to the dust Debye length, but instead takes place over the surface of a 
sphere having a much larger radius (a few times the ion Debye length). The extremely sharp 
cut-off of ip beyond f completely decouples dust grains from each other if their interparticle 
separation distance exceeds f a . 



VII. CRYSTALLIZATION 



When a > f a each dust grain is decoupled from neighboring dust grains and so the dust 
grains behave as a gas of non-interacting particles. However, if a < f a , then each dust grain 
is within the shielding cloud of its neighbor and subject to the unshielded repulsive force of 
its neighbor. Because the repulsive force scales as Zip, this repulsion becomes enormous as 
soon as a dust grain tries to move any significant distance inside of the f = f layer (i.e., 
inside of the ip — 1/Z layer). When experienced by other dust grains, a test particle dust 
grain thus acts like a hard sphere with radius f G . Thus, dust grains cannot move significantly 
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FIG. 3: Comparison of experiment parameters (short solid line marked 'expt') with model predic- 
tion that dusty plasma will crystallize if experiment parameters intersect the condensation curve 
or are above it. 



inside f Q and so the condition for strong coupling and crystallization is that a becomes less 
than f . The condensation curve is found by making the following sequence of calculations 
at each point d,fd in dusty plasma parameter space (i) calculate a(a,fd), (ii) calculate f\ 
and f Q by solving the nonlinear coupled Eqs.(60) and (61), (iii) plot the locus of the curve 
f — a and establish which side of this curve corresponds to a < f Q . The uppermost plot in 
Fig. 3 shows contours of constant f Q /a calculated for Ref. [3] and marks the contour where 
f /a = 1 as the 'condensation curve'; above the condensation curve (and with fd « 1 as 
discussed earlier) the dusty plasma should be crystallized. It is seen that portions of the 
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experiment line lie above the condensation curve which means that the model predicts that 
the dusty plasma of Ref. [3] should be crystallized. Thus, there is excellent agreement 
between the model and the experimental parameters of Ref. [3]. 

The other plots in Fig.3 are similar, but use data from the experimental results reported 
by Thomas et al. [4], Melzer et al. [6], Hayashi and Tachibana [5], and Takahashi et al. [28]. 
There is excellent agreement between the model and all these experiments with the exception 
of the Hayashi/Tachibana experiment where the experimental curve lies slightly below the 
condensation curve. The upper part of Table 1 lists the parameters of these experiments 
while the lower part gives the results of dust charging theory and then the results of this 
model using a best-fit density that is within experimental error. The main result is the 
values of fj and f . The fj values in Table I greatly exceed f d showing that the depression 
of the grain surface potential due to shielding is only a slight effect. Table 1 shows that 
when the published parameters of the experiments are used, a value of f Q is calculated which 
is slightly larger than a; the calculated ratio f Q /a is given in the bottom line of Table 1. The 
fact that f /a is greater than unity indicates that the experiment is above the condensation 
curve and so should be crystallized — this is our main result (the slight disagreement of the 
Hayashi/Tachibana experiment will be discussed later). For reference, Table 1 also lists the 
value of T associated with these experiments, and it is interesting to note that according to 
our model T has no physical significance regarding condensation and so it is not surprising 
that T has a range of quite different values for the different experiments. 

It has not been possible to compare the model to the experiment underway [8] on board 
the International Space Station, because plasma densities and temperatures have not yet 
been provided for that experiment. 

Figure 4 shows plots of logip, ip, ip on an expanded scale (to show the behavior when 
tp ~ l/Z), n e /n e0 , Ui/riio, and rid/n d0 for the Chu and I experiment [3] using the values of 
fi and f Q listed in Table 1. The tp(f) plotted in Fig. 4 is calculated using Eq.(52) in region 
1, Eq.(57) in region 2, and Eq.(35) in region 3; the electron, ion and dust densities in Fig. 4 
are calculated using Eqs.(20), (25), and (21) respectively. The dust temperature T d has been 
assumed to equal Tj so that Z = Z d ; different values of T d would only change the decay rate 
of ip outside of f , but would not change the values of fj and f since these are insensitive to 
the value of Z as long as it is large compared to unity. It is seen that there is a sharp cut-off 
of the potential at f Q and that, beyond this radius, the potential decays precipitously with 



21 



a characteristic scale length given by the dust Debye length. The potential curve is smooth 
all the way from the dust grain surface to infinity; this smoothness results from choosing fj 
and f to match ip and its derivatives at the interfaces between collisionless regions. 
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FIG. 4: Solutions for nominal parameters of Chu and I experiment as a function of r. As shown 
in Table 1, the relevant parameters are T e = 2 eV, Tj = 0.03 eV, a = 2.61, f = 0.12, fj = 1.51, 
r = 2.63, a = 0.96, and Z = 4.8 x 10 3 . From top to bottom plots are: log 10 ip, if), 10 4 ip which 
gives an expanded scale to show the region 3 decay, n e /n e o, nj/njo, and n^/n^ - The functional 
form of ip is determined from the appropriate asymptotic form in each of the three regions. Note 
the sharp cut-off of the dust density at f = f a ; the scale length of this cutoff is the dust Debye 
length \ Dd . 
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A question arises regarding why the Hayashi/Tachibana experiment [5] lies slightly be- 
low the condensation curve (i.e., has f /a = 0.97 rather than above unity). Examina- 
tion of the parametric sensitivity of the model predictions shows that intersection of the 
Hayashi/Tachibana experiment with the predicted condensation curve occurs if the assumed 
electron temperature is increased to T e = 8 eV or if the assumed grain diameter is doubled. 
Increasing the assumed atomic mass number to values larger than 16 also causes the exper- 
iment curve to approach the condensation curve but this effect is minimal in the relevant 
parameter range. Reference [5] reported an ion temperature which was not measured, but 
assumed, and an electron temperature which was estimated based on earlier measurements 
[29] made in another plasma under similar conditions. It is possible therefore that the slight 
discrepancy between the model predictions and the Hayashi/Tachibana experiment results 
from an inaccurate estimation of the electron to ion temperature ratio in the region of the 
dust grains. Better agreement would be obtained with a higher T e /T; ratio and the val- 
ues of T e and Tj used in Table I were chosen to correspond to room temperature ions and 
the electron temperature measurement given in Fig.4 of Ref. [29]. Comparison with a scan- 
ning electron microscope measurement has shown that the Mie scattering technique used 
by Hayashi/Tachibana to measure the dust grain diameter is quite accurate [30], so it is 
unlikely that the discrepancy between the model predictions and the Hayashi/Tachibana 
experiment is due to a factor of two error in measurement of the dust grain diameter. 

Another question to be addressed is the possible importance of (i) barriers due to lo- 
cal maxima in the effective potential [12, 26] and (ii) ion trapping/detrapping due to 
collisions[14-16]. The model presented here argues that the amount of net charge in re- 
gion 1 is so small that the potential in region 1 is nearly the same as the vacuum potential 
that would be produced by a bare, unshielded dust grain [see Eq.(45)]. Changing the amount 
of charge in region 1 by factors as large as order unity would not affect this argument. Thus 
any reduction in the amount of charge in region 1 because of effective potential barriers will 
make no difference to the region 1 solution, because it is already assumed that there is no 
charge in region 1. So long as the number of trapped ions in region 1 is small compared to 
Zd, the potential in region 1 is mainly due to the dust charge and again it is reasonable to use 
the vacuum potential in region 1. As for region 2, the ion density predicted by collisionless 
theory in region 2 is identical to the linearized Boltzmann relation obtained from collisional 
theory [see Eq.(37)]. Since trapping and detrapping result from collisions, trapping and 
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detrapping should tend to make the system more Boltzmann-like, but since the system is 
already Boltzmann-like in region 2, trapping and detrapping should not cause significant 
changes to the region 2 ion density profile and thus should not significantly affect the ip 
profile in region 2. Effective potential barriers in region 2 may rearrange the radial charge 
distribution in region 2 slightly, but this should cause only a small effect on ip because 
is a double integral with respect to radius of the net charge distribution [see Eq.(32)]. As 
for region 3, the normalized potential ip is so small in region 3 that ions are unaffected 
by any spatial dependence of the potential; any corrections to the region 3 potential profile 
should therefore have negligible effect on ion trajectories. Thus, while effective potential 
barriers and trapping/detrapping may modify the net charge radial profile somewhat, these 
should have a much reduced effect on the ip profile and so should not cause any substantial 
changes in the values of fj or f Q . Small changes in the ip profile should not affect the basic 
premise that there exist three concentric collisionless regions each with distinct physics nor 
the conclusion that dust grains condense when the radius f of the interface between regions 
2 and 3 exceeds a, the nominal intergrain spacing distance. 

At this point in the discussion it is possible to revisit the assumption made at the end 
of Section III that ion capture by dust grains may be ignored when characterizing the 
collisionless ion velocity distribution function in regions 1-3. Ignoring ion capture by the 
dust grain is tantamount to saying that all ions entering the collisionless region are reflected 
radially so that there are equal numbers of ions moving radially inward and outwards; if some 
ions were captured by the dust grain, there would be fewer ions moving radially outwards 
than inwards. The number of captured ions can be estimated using OML theory [31] which 
shows that the effective cross-section for ions entering from a radius where the potential 
is zero and then being captured by a dust grain is <J cap ture ~ (1 + ipdj^'d- This capture 
cross-section is to be compared to a enter = nf^ the cross-section for ions to enter region 
2 from outside (the outer boundary of region 2 is used because this denotes the edge of 
the potential well seen by the ions). Of the ions that enter region 2, the fraction that are 
captured by the dust grain is given by the ratio & capture/ Center — (1 + ipd)^d/^o- Evaluation 
of this ratio using the Chu and I parameters in Table I (ipd = 46, = 0.12, f = 2.63) gives 
& capture/ Center — 0.1 which shows that the fraction of ions captured is small enough to be 
neglected (similar results hold for the other experiments listed in Table I). This validates the 
assumption made in Section III that distinctions [20, 26] between the outward and inward 
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ion velocity distributions /+, /_ may be neglected and confirms that Eq.(17) is a suitable 
representation for the ion distribution function. 

VIII. SUMMARY 

The standard linear fluid analysis of Debye shielding fails when \q<fi/KT\ exceeds unity 
because the linear Debye shielding model is based on the assumption that \qcj)/nT\ is small 
compared to unity. This issue is important for condensation of dusty plasmas, because 
condensation requires having \Z d e<f>/ KT d \ exceed unity. 

Dusty plasmas can be characterized by an a, fd parameter space where a and fd are the 
inter-grain spacing distance and grain radius normalized to the ion Debye length. An exper- 
iment corresponds to a point in this parameter space and if the density of the experiment 
is not known precisely, then the range of densities within experimental error corresponds to 
a slanted line segment in this parameter space. 

Because shielding distances are much smaller than an ion collision mean free path, ions 
can be considered as collisionless in the shielding sphere surrounding a dust grain. A colli- 
sionless Vlasov model is used to calculate particle densities in the electrostatic potential of 
a dust grain test charge. This collisionless theory gives the same results as does Boltzmann 
theory for electrons and for dust grains because they are negatively charged but gives results 
different from Boltzmann theory for ions in the vicinity of the dust grain. Ions near a highly 
charged dust grain test particle fall into a deep potential well and are accelerated to high 
velocities. This means that no ions have zero velocity near the dust grain test particle and 
so integrals over the velocity distribution have a lower velocity limit corresponding to the 
minimum velocity of an ion falling into the deep potential well. This invalidates the fluid 
theory concept of pressure because pressure is based on the assumption of the existence of 
random velocities about some mean. For |e0//tTj| << 1, the ion density corresponds to 
the Boltzmann prediction, but for |e0//tTj| >> 1 the ion density is much less than that 
predicted by the Boltzmann relation. 

The potential in the vicinity of a dust grain test particle has three distinct types of 
behavior. These behaviors occur in concentric spherical regions consisting of (1) an inner 
region where the potential is vacuum-like, (2) a middle region where the potential includes 
both growing and decaying Yukawa-like terms with characteristic scale lengths of the order of 
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the ion Debye length, and (3) an outer region with a rapidly decaying Yukawa-type solution 
having a scale length of the order of the dust Debye length. Region 1 physics differs from 
fluid theory, is consistent with dust charging physics, and avoids the paradoxes intrinsic to 
fluid theory at large |e0//cT|. 

For any point a, in dusty plasma parameter space, the requirement for smooth matching 
of the solutions at the interfaces between the three inner regions determines the locations 
fj and f of these interfaces. Condensation occurs when a < f and occurs on the line 
a = f Q (a, f d ) which gives a 'condensation curve' in dusty plasma parameter space; a is less 
than f Q above this curve and in this region the dusty plasma is crystallized. The model 
predicts condensation parameters in good agreement with published experiments. 
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TABLE I: Comparison between model predictions and experiments 



First author & reference 




Chu [3] 


Thomas [4] Melzer [6] Takahashi [28] Hayashi [ 
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